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The present paper is concerned with the concept of the one-way quantum computer, 
beyond binary-systems, and its relation to the concept of stabiUzer quantum codes. 
This relation is exploited to analyze a particular class of quantum algorithms, called 
graph algorithms, which correspond in the binary case to the Clifford group part of a 
network and which can efficiently be implemented on a one-way quantum computer. 
These algorithms can "completely be solved" in the sense that the manipulation of 
quantum states in each step can be computed explicitly. Graph algorithms are precisely 
those which implement encoding schemes for graph codes. Starting from a given initial 
graph, which represents the underlying resource of multipartite entanglement, each step 
of the algorithm is related to a explicit transformation on the graph. 



1 Introduction 

The concept of the one-way quantum computer^ which has first been introduced by Ro- 
bert Raussendorf and Hans Briegel P El El S El E] j describes a realization of of quantum 
algorithms which goes beyond the usual network picture. A highly entangled multi-partied 
state, called cluster state P], is the basic resource for running quantum computations. The 
sequential application of von Neumann measurements to single qubits in an appropriate or- 
der, implements a quantum algorithm. It has been proven that every logical network can be 
simulated in this way jHlEI- 

For describing the preparation procedure of a cluster state, we consider finitely many 
qubits, which are arranged in a cubic lattice J = {!,■■ ■ , n}^. Furthermore, we suppose that 
an Ising interaction can be switched on and off between neighbored qubits in a controlled 
manner. For an appropriate description of this "interaction pattern" , it is convenient to 
introduce the symmetric matrix A ~ (A(i,j)|i,j S /) whose entries are A(i,j) = 1 for 
neighbored positions i,j and zero else. The preparation of the cluster state can be performed 
in two main steps. In the first step, a product state is prepared by addressing each qubit 
individually in the "standard state", which can be represented by the vector :^(|0*) + 
Here (|0'), is an orthonormal basis (computational basis) which spans the system Hilbert 
space Hi = corresponding to the qubit at position i E I . In the second step, the next- 
neighbor Ising interaction is switched on for a suitable time interval and then switched off. 
This produce, which will be explained in more detail later, creates entanglement between the 
cluster qubits and yields the desired cluster state. 

A crucial point is here, that a cluster state can completely be characterized by its stabilizer 
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group. Namely, a vector ^I'a e Hi := ®i^iHi, which represents the cluster state, fulfills a 
family of eigenvalue equations 

n<^'''^x,vi/A=±*A , (1) 

where each of the equations is labeled by a qubit position i G /. Here x.; (z^) the arc unitary 
operators on Hi, acting by Ux- (cz-) Pauli operators at position i Cz I and trivially on the 
remaining positions. The operators J^^^j z^'-"'''''xi, i G I, are called "correlation operators". 
They commute mutually and therefore generate an abelian group, the stabilizer group of 
the cluster state. The cluster state can therefore by viewed as a stabilizer code with one- 
dimensional range. 

This shows that cluster states are in fact directly related to stabilizer codes: They are 
stabilizer codes. In the usual network picture the "Clifford part" of an algorithm is composed 
of CNOT gates, local Hadamard transforms and 7r/2-phase gates. This part can be simulated 
on a one-way quantum computer by just performing measurements in eigenbases of the Pauli 
operators. An important observation is the following: Suppose one applies a local measure- 
ment in one of the eigenbases of the Pauli operators to the cluster state. Furthermore, one 
selects the state which corresponds to the "standard measurement outcome 0" for each mea- 
sured qubit. Then this state fulfills a set of eigenvalue equations, similar to from which 
follows that the received state is again a stabilizer code with one-dimcnsional range. Using the 
fact that each "stabilizer code" has a representation by a "graph code" 1101 ITT] , one knows 
already, that one obtaines a cluster state again (modulo local Hadamard transformations) 
ISj. Thus the Clifford part of an algorithm can be viewed as a sequence of transformations 
on graphs. These aspects are based on the fact that the theory of stabilizer codes is linked 
to cluster states by identifying cluster states with "graph codes". As described above, the 
symmetric matrix A ~ {A{i, j £ I) can be viewed as an "interaction pattern". But A can 
also be viewed as an "adjacency matrix" of a graph, where two vertices i,j are connected if 
A(z, j) 7^ 0. According to the equations determine the graph code associated with A. 

The theory of "stabilizer codes" is one of the most prominent techniques, for constructing 
quantum error correcting schemes |12l 1181 ICT 1151 For a review, wc rcfcrc the reader to 
\n\ . Based on the results of ^OIj the general theory of stabilizer codes from a graph code 
point of view, is discussed in Section [21 Wc provide the results which are needed for the 
analysis of cluster states and local von Neumann measurements. In particular, these results 
enable an explicit computation of the stabilizer group of the state which is obtained after 
"standard measurements" in eigenbases of the Pauli operators has been applied. 

A general description, how algorithms are implemented on the one-way quantum computer, 
is given in Section |21 Wc use here the Hcisenbcrg picture in order to stress the operational 
aspects. The fundamental resource, the cluster state, is determined by an "interaction graph" . 
Each qudit is labeled by a vertex of the graph. Those qudits which arc placed at the "input 
vertices" are prepared in the state that carries the initial quantum information. All the 
remaining qudits. positioned at the "output vertices" as well as the "measuring vertices" 
arc prepared individually in the non-binary analogue of the "standard state" . The "output 
vertices" label the positions of the qudits which carry the processed information after the 
computation has been preformed. The "measuring vertices" label the positions of those qudits 
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to which the local measurement operations are applied. The interaction graph describes one 
elementary step of the dynamics which is performed after the initial preparation has been 
completed. This is the essential operation which creates the entanglement between the cluster 
qudits: The edges corresponding to those pairs of qudits which interact. A computational 
process is now realized by performing local von Neumann measurements sequentially at the 
"measuring vertices" . 

For the implementation of an algorithm by a one-way quantum computer we have to solve 
two main tasks: 

• The "first task" is related to the fact that all measurement outcomes are completely 
random. In order to perform a deterministic algorithm, we have to "compensate the ran- 
domness" by appropriate conditional local unitary operations. This is analogous to the 
situation within a teleportation protocol. Here Bob has also to perform a local unitary 
operation conditioned on the measurement result he received from Alice. Heuristically 
the process of one-way quantum computing can be interpreted as a sequence of tele- 
portation sachems which propagate the processed information through the cluster. For 
a particular measurement result, called the "standard measurement outcome", we can 
keep the system as it is, i.e. we do not have to perform a correction by means of a local 
unitary operation. 

• The "second task" is to find a solution of the following problem: There may be different 
measurement strategies which lead to the same algorithm. Therefore, we are interested 
in computing the algorithm explicitly (as a unitary transformation or isometry) from 
the measurement strategy. Provided the problem of compensating the randomness has 
been solved, it remains to analyze how standard measurements operates on a cluster for 
a given graph. 





Fig. 1. Example for removing the measuring vertex n with self-Unk. 

A special class of algorithms, called "graph algorithms", is discussed in Section 0] These 
algorithms are based on "graph measurements" , a particular class of local measurements 
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Fig. 2. Removing two measuring vertices n, k that are connected. 



which are non-binary generahzations of measurements in eigenbases of Pauh operators. For 
these measurements, both tasks, "the compensation of randomness" as well as computing the 
algorithm from the measurement strategy can be solved explicitly. First of all, we show that 
each graph measurement can be related to an equivalent measurement strategy that only 
uses measurements in "x-basis" which is the d-dimensional analogue of the tTa;-eigenbasis. 
In view of the first task, we give an explicite construction of local unitary operations that 
compensate the randomness of measurements in x-basis. At this point we have to assume that 
the underlying interaction graph is "basic" . This property ensures that sufficient entanglement 
is created by the dynamics for processing and propagating information through the cluster. 
The second task is concerned with the computation of standard measurements in x-basis 
which can directly be computed by applying an appropriate transformation to the underlying 
interaction graph. This transformation "removes" the corresponding measuring vertices and 
builds up a new graph on the remaining vertices. For the qubit case, there are simple graphical 
roles in order to compute this transformation: 

• Suppose the measuring vertex has a selflink. Then take the subgraph by removing 
the measuring vertex. Add further edges by mutually connecting those vertices which 
are linked to the removed vertex including self links. Fig. ^ and Fig. [3 represent this 
procedure. 

• Suppose now that two measuring, without selflinks, are connected by an edge. Then we 
can proceed as follows: Build the subgraph by removing both vertices. Then add edges 
by connecting each vertex, which is linked to one of the removed vertices, with every 
vertex that is linked to the other one. (See Fig.|21for illustation). 

Unfortunately, not all sets of measuring vertices can directly be removed. This situation 
occurs, however, if one is concerned with a measuring vertex, which is not connected to another 
measuring vertex. In this case we have to apply a "local Fourier transform" (Hadamard 
transform) at a neighboring output vertex which "creates" a new measuring vertex connected 
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to the previous one. Then we can apply the corresponding removing procedure, described 
above. 

Our techniques provide some additional applications which goes beyond the analysis of 
algorithms. Examples are the simulation of interactions, quantitative analysis of multipartite 
entanglemnet, and quantum error correction: 

• Standard measurements in x-, y-, and z-bases can be represented as operations on the 
underlying interaction graph. For a given measurment strategy one obtains a new graph 
which represents a different interaction. Thus our method provides a tool for computing 
interactions which can be "simulated" . 

• The persistency of a pure multipartite quantum state is the minimal number of local 
von Neumann measurements which one has to perform to get a product state [S]. This 
is a quantity that measures the amount of multipartite entanglement. Our technique to 
remove measurement vertices can be used to compute bounds on the persistency (See 
Fig. EI). 

• Graphs, which are related to quantum error correcting codes [HI, are "basic" which 
implies that each graph code can directly be implemented by a graph algorithm. There- 
fore, by the equivalence of graph and stabilizer codes lllj . each stabilizer code can 
be realized by a graph algorithm. 




Fig. 3. The four qubit cluster state for the totally connected graph has persistency = 1. One just 
has to perform a measurement in y-basis at one of the qubits. Removing the measurement vertex 
yields a totally disconnected graph corresponding to a product state. 



2 On a general theory of stabilizer codes 

This section is dedicated to the stabilizer formalism |12lllHl[CTll()lll7j from a graph code point 
of view 201- We supply here useful techniques and results for the analysis of cluster states 
and local measurement operations performed upon. To setup the formalism, we introduce the 
concept of discrete phase space as a symplectic space over a finite field and study basic sym- 
plectic transformations. By performing the "canonical quantization procedure" , the discrete 
Weyl algebra results as an appropriate description for the observables of the quantum system. 
The symplectic transformations can be implemented by unitary operators, yielding the dis- 
crete analogue of the "metaplcctic group". The concept of stabilizer algebras (groups), which 
are particular subalgebras of the Weyl algebra, plays a central role for our analysis. They are 
naturally represented on the Hilbcrt space of the system of qudits under consideration. We 
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show here that this representation can exphcitly be decomposed into irreducible representa- 
tions, each occurring with the same multiplicity. Stabilizer codes are precisely defined by the 
corresponding multiplicity spaces. 

2.1 The discrete phase space 

The configuration space of a classical digit is an alphabet of d letters, modeled by a d- 
elementary set. As far as our analysis is concerned, we model the finite "alphabet" by the 
cyclic fields = {0, 1, • • • , d — Ijaddition and multiplication modulo d}. For binary systems 
we are concerned with the field of two elements F2 = {0, 1}. A configuration of a classical 
register is given by a set of positions I and a tuple = {q{i)\i G /) in F^, i.e. a register entry 
q{i) from the alphabet is assigned to each position i. In particular, the configuration space 
of a classical register is a liner space over the finite field F^. The corresponding phase space 
is modeled by the vector space = x F^. The first vector entry of a point {p^\q^) G is 
interpreted as momentum and the second entry as the position. The phase space is equipped 
with a canonical symplectic form. It assigns to a pair (pjlg^), of phase space vectors 

the "value" - {pi,qi) in F, with (p^,q^) = J2iei Pi^li^)- 

Symplectic (canonical) transformations are linear transformations that preserves the sym- 
plectic structure. The first type of symplectic transformations, which we consider here, is 
interpreted as "discrete dynamics" on phase space. To a symmetric matrix F — (F(i,j) = 
€ ^) associate the symplectic transformation that maps a phase space vector 
(p^k^) to {p^ ~ Tg^l?^). This "shear" transformation is viewed as one elementary step of 
evolution. The second type of symplectic transformations, that is relevant for our analysis, is 
determined by a matrix C/ — (C(i, j) G ¥d\i E I,j G J) which is invertible. That is, there is 
a matrix Cj = {C'{j,i)\i e I,j e J) such that CfCj ij and CjCj = l| holds, where l| is 
the matrix of the unit operator on F^. The matrix C/ induces a symplectic transformation 
which maps the vector {p^\q^) to {Gjq^\ — ^cjp^)? 

2.2 Discrete Weyl systems 

The Hilbert space, describing a quantum register of qudits, is the space L2(F^) of complex 
valued functions on F^ and its complex dimension is d'^' where \I\ is the number of elements 
in /. The scalar product of two functions ■01 1 V-'2 is given by 

(^1,^2) - j dq' ^i{q')Mq') (2) 

where the integration is performed with respect to the normalized Haar measure of the additive 
group F^, i.e. a function / on F^ is integrated according to 

dq' f{q') = d-l^l fil') ■ (3) 

The shift operator x{q^ ) , associated with position vector q^ , is the unitary operator, which 
translates (shifts) a function i/j G -L2(F^) by q^ . The value of the shifted function at q{ 
is {x.{q^)'ip){q{) = ipili ^ Q^)- The multiplier operator z(p^), which is associated with the 

"We denote by the transpose of a matrix Mj. 
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momentum , is the miitary multiplication operator, acting on a function ip e L2(F^) by 

(z(/)V.)(/) = xip'k'Mq') with = e^<?''-'''> . (4) 

By combining shift and multiplier operators, we associate to each point in phase space 
(p^k^) S the unitary operator 

w(/|gO=z(/)x(g^), (5) 

which we call the Weyl operator associated with {p^\q^). They satisfy a discrete version of 
the Weyl (or canonical) commutation relations 

w(p(|qf)w(p^|(7^) = xipi\qi)^ipi + pM + qi) ■ (6) 

The product of two Weyl operators is again a multiple of a Weyl operator. By taking all 
their linear combinations, we obtain an algebra 2t(/), the discrete Weyl algebra, of operators 
which act on the Hilbert space i2(F^). This algebra serves as the observable algebra of a 
quantum register. For a subset J C / of positions, the Weyl operators w{p'^\q'^) generate 
sub-algebra 2l(J) C identifying a quantum sub-register by operating only non-trivially 

on the positions J? 

The states of the quantum register, with qudit positions /, are given by the density oper- 
ators p on L2(F^), i.e. the expectation value of an operator A is evaluated by iLi{A) = tr{pA). 
One example, the so called standard state on is the pure state which is implemented by 

the projection onto the normalized vector ^[/j = ^jq/j e L2(F^). This vector is given by the 
function that assigns to each classical register configuration the constant value 1. A crucial 
prperty of the vector ^[/j is that it is invariant under all shifts x(q^) and that it is a product 
vector ^[7] = (g)ig7^[j]. By the application of multiplier operators z{p^) to ^[/j, we obtain an 
orthonormal basis, called the x -basis, 

e - { V] := z(p^)^[,] I / e F^} . (7) 

It is a joint eigenbasis for all shift operators x{q^), according to the Weyl commutation 
relations, 

x(g')C[p.] = x{p'\q'Klp'] (8) 

which immediately implies that the vector ^^i] can be viewed as an eigenvector of the mo- 
mentum operator subject to the spectral value p^ . 

Furthermore, the z-hasis is of importance here. It can be obtained by an application of 
shift operators to the normalized wave function C,[i] which is given by 

C[/](g^) = Vrf'"%^). (9) 

Since C[7] is only supported at the zero configuration q^ =0^, it is invariant under all multiplier 
operators z(p^). As for the standard state, it analogously implements a state which is now 

'The Weyl operators form a basis for the space of all linear operators on L2(F^). Hence the algebra 2l(/) is 
the full matrix algebra of all linear operators acting on L2(F^). 
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invariant under all multiplier operators z(p^). By acting with the shift operators x(g^) on 
Q/], we obtain the z-basis 

C = {C[,^] :-x(g^)C[,] |/eF^} . (10) 

Once again, it follows from the Weyl commutation relations jnj that it is a joint eigenbasis 
for all multiplier operators z(g^). 

2. 3 Implementation of symplectic transformations 

In the classical context, a symmetric matrix T implements a symplectic transformation on 
the discrete phase space which we interpret as dynamics. The quantized version of this 
dynamics is given by the unitary multiplication operator 

(u(r)V)(g'') =r(r|q^)V(q^) with T(r|g^) = e^^^'-^^^'^ . (11) 

Indeed, u(r) implements the symplectic transformation (p^q^) i~> {p^ — ^jq^lq^) i-e. the 
commutation relation 

u{T)^v{p'\q') = T{-r\q')wip' - Tq'\q')u{r) (12) 

holds. The symmetric matrix T describes a "pattern" of two-qudit interactions. This can be 
visualized by a weighted graph whose vertices are the positions i G /. Two vertices i,j are 
connected by an edge if the matrix element T{i,j) 7^ is non- vanishing. The value of the 
matrix element T(i,j) is interpreted as the "strength" of the interaction. 

The symplectic transformation {p^\q^) i— > {Cjq^\ — *Cjp^) is also implemented by a unitary 
transformation F^^cj] ■ called the Fourier transform associated with the invertible matrix 
C/. It identifies the Hilbert space L2(F^) with L2{¥'^) according to 

{Fic,']^){q') = J ,Wx{Ciq'\q'mq') . (13) 

By construction, the commutation relation 

F[c/]w(/|g^) = w(C/g^| - *C,V)^^[c/] (14) 

follows. For our purpose, we consider those Fourier transforms which are associated with an 
invertible matrix C/, where / and J are disjoint. We also visualize such a Fourier transfor- 
mation by a weighted graph on the union / Ji'^A vertex i € / is connected with a vertex j € J 
if C{j, i) 7^ 0. The symmetric matrix C/ + *Cj is the adjacency matrix of the graph. If each 
vertex of / is connected with precisely one vertex in J, then we call C/ a connecting matrix 
or connecting graph to put emphases on the graph-theoretical aspects. The corresponding 
Fourier transforms play a crucial role. They are local unitary transformations on the qudit 
systems labeled by J. In fact the Fourier transform, associated with a connecting graph, is a 
tensor product 

^[c/l=(8)^^[cPi (15) 
je/ 

'^Notation and conventions; For the union of two finite sets K and L we often drop the union-symboh KL := 
KUL. 
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of Fourier transforms operating independently on each single digit. Namely, a connecting 
graph Cj decomposes into connected components 

(16) 

where is the connecting graph consisting of the line with weight C(j, i) connecting the 
vertex i €z I with its unique counterpart in J. 

2.4 Stabilizers 

An isotropic subspace of a discrete phase space is characterized by the vanishing of the sym- 
plectic form for all pairs of its vectors. Making use of the results of JHIi an isotropic space can 
be parameterized by a weighted graph A on a union of four sets the input vertices /, output 
vertices J, auxiliary vertices K, and syndrome vertices L. The isotropic space, associated 
with A, is the subspace which is defined according to"* 

S[ij,K\A] {iA',K<l'''W') I ^jW = 0} . (17) 

According to the Weyl commutation relations the set of Weyl operators 

{wiA''jj,q-"'\q'') I q''^ E kerA^^} . (18) 

generate an abelian sub-algebra 2l(/, J, A'| A) C 2l( J) which we call the stabilizer algebra of the 
isotropic subspace j_^|a]. The stabilizer group is just the unitary group generated by H18|) . 
The definition of the isotropic space as well as the construction of the stabilizer algebra only 
depends on the sub-graph Ajj^ subject to the complement of the syndrome vertices UK. 
Despite of this, the edges that connect syndrome vertices with others play an important role 
for the representation theory of the stabilizer algebra. 

However (see \W\). for parameterizing an arbitrary isotropic space we can restrict ourselves 
to weighted graphs which we call admissible. They fulfills the following list of conditions: 

(Gl) The set of vertices satisfy |/| + \L\ = |J|. 
(G2) The submatrix Ajj^j^ is invertible. 

(G3) There are no edges that connect syndrome vertices, i.e. A£ = 0. 

Simple examples for admissible graphs are shown in Fig. 01 As shown in ^ . these graphs 
correspond to [[5,1,3]]^ stabilizer codes. However, one only has to consider the syndrome 
vertices if one is interested in an explicit construction of the decoding operation. The encoding 
itself just makes use of the subgraph which is obtained by wiping out the syndrome vertices. 

Since 2l(/, J, K\A) is an abelian C*-algebra its representation on ^2(1^^), which we call its 
natural representation, can be decomposed into irreducible representations, called characters^ 
Each character occurs with a certain multiplicity, i.e. the Hilbert space L2{¥'^) decomposes 

''Notation and conventions: For a matrix S^j and for two subsets K C M, L C N, wc write Sj^ = k)\l G 

L,k £ K) for the corresponding sub-block. 

'^An irreducible representation of an abelian C*-algebra C is one-dimensional, i.e. it can also be viewed as 
a state ro : 1— » taj{A) on C The functional zu is multiplicative and therefore a pure state. 
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Fig. 4. Examples for graphs, corresponding to [[5,l,3]]ti stabilizer codes, fulfilling the conditions 
(G1)-{G3). 



into a direct sum of multiplicity spaces. The decomposition of the natural representation into 
irreducibles can directly be derived from the structure of the weighted graph A. In partic- 
ular the syndrome vertices are used to arrange the combinatorics for obtaining an explicit 
decomposition. How this can be done, will be discussed in the next section in detail. 

2.5 Representation theory 

The stabilizer algebra 2t(/, J, K\A) is generated by the set of Weyl operators H18|l . They form 
a linearly independent family of operators and since the product of two of them is a multiple 
of another Weyl operator within this family, it is a basis for J, ^| A). The condition (G2) 
implies that the kernel of Aj^ contains d'^' elements and the dimension of 2t(/, J, isTjA) is 
therefore d'^L It turns out that for all characters of the stabilizer algebra the multiplicity 
spaces have the same dimension. And since the dimension of i2(IP'd) follows from (Gl) 

that the each multiplicity space is d'^' dimensional and therefore isomorphic to L2{¥^^). There 
are d'^' different characters of the stabilizer algebra and we may label each multiplicity space 
by a classical register configuration in F^. The basic idea to perform the decomposition is 
to construct for each classical register configuration an isometry from L2{¥^^) into ^2(1'';^) 
whose range is precisely the multiplicity space labeled by q^. This is nothing else but a coding 
operation. The Hilbert space L2{¥^^) is the input space corresponding to the system we wish 
to encode, and L2{¥'^) is the output space in which quantum information is encoded. This 
point of view is the reason for calling the set / input vertices and the set J output vertices. 

The x-basis vectors ■^[pW] and its z-basis counterparts C[q'^] = ll^) induce isometric map- 
pings 

:7^^L2(F5')®H, (19) 
n[^«] :H-*L2(F^)®H, i^^C[g«]<»V' (20) 

were we write ^[n] = *i'[o"] ^-nd IIj^v] = n[o™]- From these isometrics and the unitary u{A) 
which implements the dynamics for the "interaction pattern" A, we build the linear operator 

V[A|,i,] := %/d'"'^'$f^^^]w(A)$[jK]n[,i,] (21) 
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that serves as a candidate for identifying the input space L2(F^) with a muhiphcity space. 
One also obtains a more explicit expression for the operator 'V[A\qi-] by applying it to a vector 
Clq'] = 1 9^) of the z-basis for the input space and expanding the result in terms of the z-basis 
of the output space 

V[A|,.]|g^) = Vd^'""^ J d/^ r(A|g^^^^)|/) . (22) 

Let us suppose that the graph A has no edges that connects a vertex with itself, then the 
phases r(A|g-^-^^^) can be obtained by just taking the product over all edges of the graph 
A, where to each edge the phase exp(^g(i)A(i, j)g(j)) is assigned. Comparing the 

expression H22(l with the construction of graph-codes, outlined in [HIEIIj '^^ see that V[A|gi] is 
nothing else but the encoding isometry of a graph-code whose range is indeed a multiplicity 
space of a character of the stabilizer algebra 2l(/, J, K\A). 

The following central theorem states that these intuitions are correct and it shows how the 
decomposition of the natural representation of a stabilizer algebra is encoded by its associated 
weighted graph. 

Theorem 1 Let A be an admissible graph on the set of vertices IJKL. Then the natural 
representation of the stabilizer algebra 2l(/, J, K\K) decomposes into a direct sum of characters 
{ra7[gt]|(7^ G F^}. The character VJ\^qL^ is uniquely determined by the Weyl expectation values 

tn[,.](w(Ai^g^^|g'^)) = r(A|g^^^) , with q'^ € kerA^f . (23) 

The operator V[yY|qi,] is an isometry whose range is the multiplicity space of the character 

Proof. A proof of the theorem can mainly be derived by using the techniques which 
has been developed in pi IIUL [THj . To keep the paper self contained, we prove some useful 
technical lemmas in the appendix. First of all. Lemma 1 121 ensures that the operators v^^igi] 
are isometric. Lemma lT^ states that the range of the isometry 'V[A\qi-] is in fact the multiplicity 
space of the character W[qL-^ of the stabilizer algebra 2l(/, J, i^|A). To conclude the prove of 
the theorem, we apply Lemma [T^ from which follows that the ranges of the isometrics v^a^^l], 
q^ G F^, are mutually orthogonal and that they "add up" (as a direct sum) to the full Hilbert 
space L2{¥j) □. 

As an abelian algebra of dimension d'^', 2l(/, J, i^|A) is isomorphic to the abelian alge- 
bra of complex function on the "syndrome configurations" . Due to Theorem ^ we can 
explicitly identify operators in 21(7, J, i^TjA) with functions on F^. Namely, a Weyl operator 
w{Ajj^q'-^^\q'^), q^^ G kerAj^, is just identified with the function q^ ^ T{A\q^^^). Fur- 
thermore, Theorem n implies that every vector ■(/; in the range of the isometry ^[A\q^\ fulfills 
the eigenvalue equation 

w(A^^g'^^|/)V^==r(A|/^^)V' (24) 

whenever the classical register configuration q^^ satisfies the constraint Aj^q"^^ — 0. This 
is just what corresponds to the eigenvalue equation for a qubit cluster state and there are 
two special situations of particular interest: 
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No input vertices: The first special case is concerned with admissible graphs for which the 
set of input vertices is empty / = 0. In this case, the vector 

*[A|,L] Vd^'"^^'t>li^^u{A){^[jK] ® C[g^]) (25) 

is the unique normalized vector in the one-dimensional range of V[A|gi.]- The stabilizer algebra 
2t(J, isr|A) = 2t(0, J, -ftr|A) is generated by Weyl operators which correspond to the isotropic 
subspace of phase space vectors {Ajj^q'''^\q'^), q-^^ e kerA^^. An immediate consequence of 
Theorem n is then the corollary: 

Corollary 2 The family of vectors {^'[A|g^-] G '^'^ orthonormal basis in L2{¥'^) and 

each basis vector ^^j^\gL^ fulfills the eigenvalue equations 

w(A:;^/^^|q^)vI/[A|,.] - r(A|g'^^")vI/[A|,.] (26) 

for all classical register configurations q'^^ E kerAj"^. 

No input and auxiliary vertices: The second situation, being concerned here, is even more 
special. Namely, we now consider admissible graphs with no input and auxiliary vertices 
I = K = %. Since there is no edge that connects a vertex in JL with a vertex in 0, the 
kernel of A^^" = is just the full vector space The corresponding stabilizer algebra 
2t(J|A) := 2t(0, J, 0|A) is generated by the Weyl operators •w{Njq'-'\q-^) where q'^ can be any 
classical register configuration. Hence, the vectors of the orthonormal basis {^'[A|qi]k^ G F^} 
fulfill the eigenvalue equations 

^{KWW')-^[A\q-] = r(A|/^)vI/[A|,.] (27) 

for all q'^ . 

3 One-way quantum computing 

The results on stabilizer codes, discussed in the previous section, are now used as a basis for 
investigating the concept of the one-way quantum computer. The algorithms under consider- 
ation are given in terms of elementary operations which are preparation of individual qudits, 
performing one elementary steps of dynamics, applying local von Neumann measurements 
and performing local unitary operations, conditioned on measurement results. 

3.1 Elementary operations 

For the implementation of algorithms, we discuss here those operations which are regarded 
as elementary. Namely, local preparations of individual qudits, elemetary steps of a global 
dynamics, local measurements and conditional local unitary operations. 

Local preparation of individual qudits: The input qudits at position / are prepared in the 
input state, carrying the quantum information one wishes to process. The qudits at the 
remaining positions N are individually prepared in the standard state. In the Heisenberg 
picture, this operation is described by an encoding channel i?/, which maps an observable 
A e 2l(/iV) to the operator 



(28) 
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If the set / = is empty, then the operation Eijj just prepares the standard state which is 
implemented by the vector C[n]- To point out the analogy to the binary case we expand it in 
terms of the z-basis, where we sometimes write the basis vector C[g"] ^ "kct" \q'^): 




For (i = 2, we indeed obtain a tensor product of the vectors -^(jO") + |1")) over N. 

Elementary step of discrete dynamics: One elementary step of the discrete dynamics, is 
determined by a weighted graph F on a set of vertices N. It is given by an automorphism 

ar : 2t(iV) ^ 2l(iV) (30) 

that maps an operator A £ 2t(A^) to 

ariA) = u{T)* Au(r) . (31) 

The dynamics operates parallel on all qudits in the cluster within one step and is, in com- 
parison to the usual network picture, a very fast operation. As already mentioned, the graph 
F describes the type of two qudit interactions: Two qudits interact if their positions are 
connected by an edge. In view of this F is called an "interaction graph" . 



Local von Neumann measurements: A von Neumann measurement, acting on the qudits 
at positions M, is determined by a choice of a orthonormal basis [3 ~ {i3[qM]\q^^ £ ¥^'} 
of L2{¥^^). In the Heisenberg picture, this operation maps a function / of the classical 
observable algebra £(Af ) of functions on F^-'^ to the operator 

%M](/):= E e2l(M) (32) 

where P[p\giii] is the projection onto the basis vector (3[gM]. Von Neumann measurements are 
mathematically characterized by the property that they preserve the multiplicative structure, 
in other words, P ~ -P[/5.a/] is an C*-algebra homomorphism. This means that the product 
fg of two functions is mapped to the product of operators P{f)P{g) and the complex conju- 
gate function / is mapped to P{f) = P{f)* to the adjoint of P{f)- A quantum state ui on 
2l(M) is mapped to the classical probability distribution. The probability for receiving the 
measurement outcome q^^ is then given by the expectation value uj{P[p^gM-^). The values for 
the classical measurement outcomes under consideration are classical register configurations. 
In the binary situation one receives, for example, a bit string of length \M\. The additive 
structure of register configurations, allows to choose the zero-register configuration q^^ = 0^^ 
as a natural "standard" for the "initialization" of a classical register. In view of this, we call 
q^^ = 0^^ the standard measurement outcome and the corresponding selective measurement 
operation is called standard measurement operation. A local von Neumann measurement is 
related to a product basis /?, i.e. each vector of the basis is a product vector. The correspond- 
ing local measurement operation decomposes into a tensor product P[f3,M] ~ ®meJ\/^'[/3,m] j 
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addressing each qudit individually. The x-basis ^ = {^[p'^'']\p^^ G F^^} as well as the z-basis 
( = {C[qA/]|(7^^ G F*-'^} consists of product vectors. Thus they correspond to local von Neu- 
mann measurements P[^,m] Sind -P[f ,m] i respectively. Since the a;-basis consists of eigenvectors 
for all shifts, the measurements in x-basis are invariant under shift operations, i.e. 

x(g'0%M](/)x(g*^)*=%M](/) (33) 

is valid for all q*^ and for all functions / e €{M). Analogously, the measurements in z-basis 
are invariant under all multiplier operations. Measurements in x- and z-bases are mostly 
applied to a subsystem 2t(A/) C 2t(/M). In this case, the measurement operation Z*!^ j\f j 
only operate on the positions M and leave the rest of the system alone. Keeping in mind 
that for all operators A E 2l(/) the identity $[pA/]^$j'pjjj = A® P[j|pA/] holds, it is a straight 
forward observation that the co-isometries {<I>*pM] Ip^''^ G F^^} are Kraus operators for the 
local measurement operation P[^,m] i where 'J'J'pM] corresponds to the random measurement 
outcome p^^ . Analogously, if one is concerned with a local measurement operation in z-basis, 
then the Kraus operators are given by the co-isometries {HJ'^jjj jg^-'^ e F*-'^}. 

Conditional local unitary operations: Let M be a set of positions at which local von Neumann 
measurements have been applied and let J be the set of positions of the remaining qudits. A 
conditional local unitary operation is given by a channel 

C : 2l( J) ^ €{M) ® 2l( J) (34) 

which is implemented by a family of local unitary operators {s[(7|pA/] in 2t(J) according 
to 

[C(A)](/0-sfc,|pM]As[c|pM] . (35) 

Indeed, this is nothing else but performing the applying local unitary S[p|pA/] in case the 
measurement result p^^ was received. 

3.2 Algorithms 

An algorithm is implemented by applying a sequence of elementary operations, as described 
above, in an appropriate order. Each algorithm is mathematically modeled by a channel, i.e. 
completely positive unit-preserving map, 

T : a( J) ^ 2l(/) (36) 

that maps an observable of the output system 2t(J) to an observable of the input system 2l(/), 
i.e. concerning the Schrodinger picture, a state of the input system is mapped to a state of 
the output system. Using the concept of the one-way quantum computer, the channel T is 
realized by the following procedures: 

Step )00): In the first place, the qudits with positions at the input vertices I are prepared in 
the state that carries the quantum information wc which to process. The qudits at the output 
vertices J and the measuring vertices K are individually prepared in the standard state. This 
elementary operation is given by the channel 

Ei : %{IJK) 2l(J) . 
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Step )01): The next elementary operation is an application of one elementary step of the 
discrete dynamics 

ar ■■ WJK) %{IJK) 

associated with an interaction graph F with input vertices /, output vertices J and measuring 
vertices K . 

Step )10): A local von Neumann measurement is applied to the input vertices / in x-basis 
and to a subset Ki <Z K measuring vertices in the product basis /3i of L2 (F^^ ) : 

PmPiJK^] ■■ €iIK,) ® 2l( JA'i) ^ 2l(/JA') 

where = K \ Ki is the complement of Ki in K . This measurement produces a random 
outcome p^^^ . 

Step A conditional local unitary operation, subject to the outcomes of the previous 

measurement, is performed: 

Ci : ■^{JK^) €{IKi) ® %{JK^) . 

It is implemented by the family of local unitary operators {s[Ci\p'Ki]\p^^^ G F^^^} which are 
localized in JK^, i.e. they belong to 2t(Ji4'^). 

Step )tO): In the tth measurement round a local von Neumann measurement operation with 
respect to a product basis f3t of L2(FJ') is applied at the measuring vertices Kt C K*~^: 

Here AT* = K*-~^ \ Kt is the complement of Kt in K*^-^. 

Step )tl): A local unitary operation, conditioned on the outcomes of the previous measure- 
ment, is performed: 

Ct : %{JK*) ^ <L{Kt) (E) '^{JK^) . 

The wanted algorithm, is realized by proceeding the scheme inductively, until all qudits, 
placed at measurement vertices, are removed. Thus the channel T is given by 

T = EioaroFoo---oFn , (37) 

where Ft := Pip^^Kt] ° is the composition of two elementary operations, a measurement and 
a corresponding conditional local unitary operation. 

The composition )00)01) of the initial steps )00) and )01) is interpreted as the "reading in" 
of quantum information into the cluster. This information is "spread out" through all qudits 
via the entanglement that is created by the dynamics. In fact, if the graph F is connected, 
then the received states is "highly entangled" in the sense of the persistence of entanglement 
[3] . This entanglement can be used to run the quantum algorithm by an sequential application 
of local measurement operation. 

The composition of the steps )00)01)10) (also symbolized by )00 10)) yields a channel 
which has an implementation by the Kraus operators 

{ ^lu<,^uMm[JK] I p'""' 6 F"^i } , (38) 
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where ui is a local unitary operator that transforms the z-basis into the basis /3i. This channel 
is not a pure quantum channel since it produces classical measurement results p^^^ which 
are completely random. In order to "control the procedure" of information processing, we 
have to require additional constraints on the conditional local unitary operation Ci which is 
applied by step Controlling a procedure means here that we have no loss of information 
during the process. This requires that the channel T is pure. By construction, the channel 
T)oo— ►11) which is given by the composition of the steps )00 11) is implemented by the 
Kraus operators 

{ S[c,ipiK,]'Py,,,^u,uim[JK] I p"''' e F^^^'i } . (39) 

Purity implies that all elements in this set are multiples of one isometry. This is the case 
if the action of Ci "compensates the randomness" the effect of the measurement operation 
Pl^^0i,iKi]- Whenever the measurement result p^^^ is received, the conditional unitary oper- 
ations rotates the system back into the state which corresponds to the standard measurement 
outcome p^^^ = 0^^^. In other words, the implementation of T^oo^ii) is "gauged" with re- 
spect to standard measurement results. Expressed in terms of the Kraus operators (|38|l and 
(|39|l this just means that the family of local unitary operators {s[Ci\p'Ki]\p^^^ G F^^^} has 
to fulfill the identity 

^[C^\p"<l]%IK,]Ulu{T)<^>[JK] = ^ljK,]Ulu{T)<i>iJK] (40) 

for each measurement outcome p^^^. As a consequence the channel T^qo^ii) can be repre- 
sented by the single isometric Kraus operator 

^/rf'"^''*[/A',]^il"(^)1>[JK] • (41) 
In order to realize the algorithm T, we are faced with the following two main tasks: 

Task 1: For each measurement round t we have to compute the local unitary operators sjc'^].] 
that compensates the randomness of the corresponding measurement. 

Task 2: The Kraus operators (|41|l . which are related to standard measurement results, have 
to be computed. 

Remark: The composed channel Fi := P[ui^jKi] ° C'l can be interpreted as adjusting the 
basis for the next measurement round -P[/32,-fs'2] subject to the received measurement result of 
the previous round. More precisely, depending on the measurement result p^^^ of the first 
measurement round, the following local von Neumann measurement operation is performed 
in the adjusted product basis S[c'^|p/Ki]/32. 

If we view for the tth measurement round ^{IKt) as Alice's observable algebra and %{JK*) 
as Bob's, then the channel T^oo— >ti> can be viewed as a teleportation scheme which propagates 
information through the cluster. After all qudits at K has been measured, the information has 
been processed via "teleporting" it from the input qudits into the output qudits. The entan- 
glement within the cluster has been used up after the computation has been performed. This 
is the reason why the computational model, described above, is called a "one-way quantum 
computer" . 
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4 Graph algorithms 

In the previous section the general scheme for implementing an algorithm via one-way quan- 
tum computing is discussed. The basic initial procedure is the initial preparation followed by 
an elementary step of the dynamics )00)01) which creates the basic entanglement resource. 
In principle the measurement operations, that follow, can be performed in any basis. The 
present section restricts to a particular class of bases that arise from stabilizer codes. They 
can be represented by graphs and the corresponding von Neumann measurements are called 
graph measurements as it is described in detail in the first part of this section. 

The algorithms under consideration make use of one local graph measurement round as it is 
described by Step )10). The first task is to find an explicit construction for a conditional local 
unitary operation that compensates the randomness of the graph measurement. It turns out 
that graph measurements can be converted into measurements in x-basis by an appropriate 
change of the underlying interaction graph. This enables for a particular class of graphs, 
which we call "basic", to solve the task. In view of Step )11), a family of Weyl operators, 
implementing the desired conditional local unitary operation, can be constructed explicitly. As 
a direct consequence, each encoding procedure of a graph code (and therefore each stabilizer 
code) can be implemented by a one-way quantum computer ^ I1U| . 

Moreover, the second task can be solved, which consists in computing the Kraus operators 
for standard measurement outcomes. We show that standard measurement operations in 
x-basis can be viewed as operations on graphs. This leads not only to the fact that graph 
algorithms can be represented by one graph, which only has input and output vertices, we 
also obtain an explicit construction scheme that computs the final interaction graph from the 
underlying data, i.e. the given interaction and measuring graph. 

4-1 Local graph measurements 

A particular class of graph measurement operations ar the so called graph measurements. 
They are related to bases which emerge from stabilizer codes via Corollary |2 in the following 
manner: We consider an admissible graph A on the disjoint union MKL of output vertices 
M, auxiliary vertices K and syndrome vertices L. Due to Corollary |21 the set of vectors 
{\E'[A|qt] Ig'^ € F^} is an orthonormal basis of the Hilbert space L2{V^-'). The von Neumann 
measurement operation Pa, operating on the system 2t(M), corresponds to this basis. It is 
the completely positive map which assigns to a classical observable / £ £(-^) the operator 

PAif) E /('^'') (42) 

in 2t(Af), where P[A|qf-] is the projection onto ^'[Aiqi]- We also get a Kraus representation 
of Pa by the co-isometries {'^[Aiq^-j 1?^ ^ where the operator is related to the 

measurement outcome q^. Note that for an observable A G 2l(/) of a system, that is not 
affected by the measurement operation, we have '^[A\q^]^'^*A\qL] = Ai^ P[A|gi]. 

Since the sets M and L contain the same number of elements, the set MKL admits a 
partition into mutually disjoint sets 



MKL = IJ {to} U k{m) U {^(to)} 

me A/ 



(43) 
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where m i~-> l{m) is a bijective map which identifies the output vertices with the syndrome 
vertices and fc(m) is a set of auxihary vertices. Suppose now the weighted graph A is a sum 
of connected components A ~ X^Mem "^here each A^ is the admissible sub-graph on the 
vertices {m} U k{m) U {l{in)}. Then the vector decomposes into a tensor product 

*[A|g^] = *[A,„|g'(™)] (44) 
me A/ 

and Pa is a local von Neumann measurement operation. The graph A,„ determines the 
measurement basis at the position m. 



Syndrom vertices O 



Auxilliary vertices 



6 



Output vertices 



Fig. 5. From left to right: The measurement graphs for local von Neumann measurements in x-, 
J/- typo- and 2-basis. 

For analyzing graph measurements, it is sufficient to study graphs on two-elementary set 
{to, 1} (one output and one syndrome vertex) or three-elementary sets {to, k, I}. For the case 
that there are no auxiliary vertices, the simplest example for a measurement graph is given 
by the adjacency matrix 

(45) 



X = 





1' 







which represents the graph consisting of one single line that connects the vertex j with I (see 
Fig. EJ. Here IJ^ is the 1 x 1-matrix block which corresponds to one edge that connects the 
vertices n and to. Obviously, X is admissible. The orthonormal basis, that consists of the 
vectors ^[xiq'j — 9' £ I*", is just the x-basis since 

*[xi,'](9'") = x(irVk'") = e[ir,'i('7'")- 

The corresponding stabilizer algebra ^{m\X) is generated by shift operators w(Ar™(7™|(j™) 
x(g'"). 

A further example, which is almost as simple is given by the adjacency matrix 



(46) 



Y = 



n 
1' 



1 



(47) 



of the admissible "tadpole" graph consisting of one line that connects the vertex k with I and 
n self-links at the vertex to (see also Fig. EJ. In this case, the orthonormal basis, which is 
given by the vectors ^'[y|gi], € F, is a y-type basis: 



*m,'](9'")=r(n|g'")x(irVir). 



(48) 
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The corresponding stabilizer algebra 2l(m|y) is generated by Weyl operators w{nq"'-\q"'-), i.e. 
the vectors 4'[y|gi], e F, form a basis of eigenvectors for the Weyl operators w(ng™|g™). 
Note, that for 77, = we obtain the a;-basis. 

The following example makes use of an auxiliary vertex. The adjacency matrix 

/ IT \ 
Z = It If (49) 

V 4 ; 

corresponds to the measurement graph which connects m with k and k with I (see Fig. EJ. 
Clearly, Z is admissible. Now, the orthonornial basis, consisting of the vectors ^'[z|g'] = 
C[i;"g']i 9' '= is the z-basis since 

vi/[^i,,](g") ^Vdjdq' xiq'l^ir + If?') = C[ir,'](9™) ■ (50) 

Indeed; the kernel of the block 

ZtmM = ( 1™ ) (51) 

consists of those vectors g^'"-*^!, for which the component g'" ~ vanishes. Thus the corre- 
sponding stabilizer algebra 2l(m, k\Z) is generated by shift operators w(l™(7*''|0™) = z(l™(7'''). 
As a consequence, all measurement procedures in a;-, j/-type and z-bases are covered by con- 
sidering graph measurements. 

Now we consider an applicaion of a graph measurement for the first round (Step )10)). 
The corresponding algorithm is called a graph algorithm. It is constructed from a interaction 
graph F, with input vertices /, output vertices J and measuring vertices M, and from a local 
measuring graph A which lives on M, the auxiliary vertices K and the syndrome vertices L. 
The channel T^^ia], which describes the algorithm, is defined by the following sequence of 
elementary operations (See Fig.E)): 




Fig. 6. The graph algorithm for an interaction graph F and a measurement graph A is represented 
on the left hand side. The sequence of the single elementary steps )00 — » 11) is indecated below 
the figure. Both four step procedures )00 11) and )00 11 ) yield the same operation. It 
realizes the graph code associated with the graph F — Z^A. 



Step )00): The quantum information, we wish to process is encoded via the operation Ej, 
i.e. each qudit at the positions IM is prepared in the standard state. 
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Step )01): An elementary step of the dynamics ar is applied, where the graph F lives on the 
union of input vertices /, measuring vertices M, and output vertices J. 

Step )10): The input qudits are measured in x-basis P[^j]- At the measurement vertices M, 
a graph measurement Pa is performed for a measuring graph A, with measuring vertices M, 
auxiliary vertices K and syndrome vertices L. This produces measurement outcomes ,q^. 

Step A conditional local unitary operation C[r,A] is applied that compensates the 

randomness of measurement P[^.i] (E) Pa- 

Let us first concentrate on the channel 7)00)01)10) which is given by composing the steps 
)00)01)10) (See Fig. IHlfor illustration). It has a representation by the set 

{ Vd"'"v[A|,.]U[r|p.] |g^eF^} (52) 
of Kraus operators. For convenience, we have introduced here the encoding operator 

V7i '$fp„]«(F)$[jM] , (53) 

which is associated with the graph F and the measurement result p^ , where A'^ is a subset in 
IM that consists of input and measuring vertices. As we shall see later, the encoding operator 
associated with F for the standard measurement outcome U[r|7v] •= U[r|o"] plays a particular 
role: The encoding operator, associated with the graph F and the standard measurement 
outcome, is the Kraus operator that implements the quantum code associated with F in the 
sense of pnT)]. 

So far everything can be expressed explicitly in terms of the given interaction graph F 
and measuring graph A. The problem, which occurs next, is to given an explicit construction 
for the operation Cp.A]- At least we have to show the existence of local unitary operators 
Sjp A|p^,<ji] which compensate the measurement outcomes p^ ,q^: 

V[A|o^]"[r|/] = ^[r.A\p' .qi']^[A\q^]^[r\p'] ■ (54) 

If we can solve the Equation l|54|l for each measurement outcome p^ , then we know that 
the channel PfriA] has a representation in terms of the single isometric Kraus operator 

Vd y^lA\o^]U[r\i] . (55) 

In order to construct explicit solutions of the equations H54|l , we make use of the fact that 
every graph measurement can be realized by just performing measurements in x-basis. More 
precisely: 

Theorem 3 For each measurement outcome {p^,q^) the identity 

|L| 

^[A\qi']^ir\p'] = V d U[r-ZLA\pi-Al"^q^] (56) 

holds, where Z^A = Aj|^ is the sub-graph which is obtained via removing the vertices L. 
Proof. Lemma provides a useful relation between the isometry V[A|<ji.] with the encoding 
operator U[p|p7]. Namely, 



V[A|gi]"[r|p^] = ^[qL]Ulr-A\0''"<,p'] 



(57) 
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holds, where F — A is a graph with input vertices /, measuring vertices KM and output 
vertices JL. Due to Lemma [T71 it follows that encoding operator, associated with the graph 
r — A and the measurement outcome (0^^^ ,p^) fulfills the identity 

U[r-A\o''"<y] = 'r(A|(7^)w(A|;q^|q^)u[r_A|p/_A£"^q^-] • (58) 

We apply the co-isometry n*qt] = nj'^]x(— g^) on both sides of this equation and we exploit 
the fact that there are no edges of A which connect syndrome vertices L. This yields the 
relation 

"^[A|gi]"[r|p^] = ^[L]^lr-A\p' -A'^"' qi^] ■ (59) 

One observes that for each function i' € L2(F^-') the identity (EJ'^j ■(/;)('7'') = Vd ^^^i'{0^ ,q-') 
holds which finally implies that 

^[L]^lr-A\p' -A'^^' qi-] = ^ U[r-ZtA|p^-Af ^^gi-] (60) 

is true. This concludes the proof of the theorem □. 

One immediate consequence of the theorem is that the channel T[r|A] can be decomposed 
by an alternative sequence of elementary operations which only include local measurements 
in x-basis (See Fig. El: 

Step )00): The quantum information, we wish to process is encoded via the operation Ei, 
i.e. each qudit at the positions JKM is prepared in the standard state. 

Step )01): An elementary step of the dynamics ch^y-ZlA] is applied, associated with the graph 
F — Z^A on the union of input vertices /, measuring vertices MK , and output vertices J. 

Step )10): A measurement in x-basis P[^jmk] is performed at the input vertices / and the 
measuring vertices MK, producing a measurement outcome p^^^^ . 

Step The randomness of the previous measurement is compensated by a conditional 

local unitary operation Cp-ZtA]- 

To verify the equivalence of both procedures, we make use of the observation, that the 
composition of the steps )00)01)10) has a representation by the Kraus operators 

{ Vrf""'''"u[r_z.A|p--] I /^^^ e F^^^^-} (61) 

each corresponding to the measurement outcome p^'^^^ £ F^^^^. If we assert here that 
the randomness of the measurement outcomes can be compensated, then for the channel 
T[p_2^A]j which is given by composing the steps )00 —^ is represented by the single 
isometric Kraus operator up-ZtAl/A/K]- By Theorem 13 this isometry coincides with the 
Kraus operator (|55|) and therefore both procedures yields the same result (Fig. right hand 
side): Tp-z^A] = 2^[r|A]- Therefore it is just one graph F — Z^A needed for determining a 
graph algorithm. 

4-2 Compensating the randomness 

As discussed in the previous subsection, all algorithms which are based on graph measure- 
ments can be converted into algorithms which arc only faced with measurements in x-basis. 
Without loss of generality, each graph algorithm can be represented by one graph F with 
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input vertices /, measuring vertices K and output vertices J. The corresponding channel 
Tr is then constructed by composing the four elementary operations )00 11) given in the 
previous section, where the graph F — Z^A is just replaced by F. 

Since we are concerned with x-basis measurements we expect, following the analysis of 
that the local unitary "byproduct" operators, which implement Cr, are Weyl operators. In 
view of this we introduce the following notion: A linear map which assigns to each mea- 
surement outcome p^^ a phase space vector (translation) Qp^^ € S'' is called compensating 
for the measurement operation in x-basis at IK if the equation 

niniK] = A(p^^)w(ep^^)*U[r|p.^] (62) 

holds for all measurement outcomes p^^ , where X{p^^) is some phase depending on p^^ . 

If we assume that there exists a compensating linear map Q, then the channel Cr is 
implemented by the Weyl operators 

S[r|p.K-] = w(ep"^) . (63) 

But what can one say about the existence of compensating linear maps yielding appropriate 
conditional local unitary operations? As already discussed above, it follows from the fact 
that the linear map Q is compensating, that is implemented by a single isometric Kraus 
operator, namely the encoding operator U[r|/K]- Thus the following proposition is a necessary 
condition: 

Proposition 4 // there exists a compensating linear map for the measurement in x-basis at 
IK , then the encoding operator U[p|/;f] is isometric. 

Sufficient conditions for the existence as well as an explicit construction of compensating 
linear maps, can be formulated in terms of the interaction graph: An graph on the union of 
input vertices /, measurement vertices K, and output vertices J, is called a basic graph if the 
block matrix Tf^ is injective. The set of all basic graphs F is denoted by BG(/, J; K). 

We point out here that the condition for a graph F to be basic implies that the encoding 
operator U[p|/;^;] is an isometry and vice versa. We show that the following is true: 
Theorem 5 Let F be weighted graph on the union of input I , output J , and measuring 
vertices K . There exists a compensating linear map Q for the measurement in x-basis at IK 
if and only ifV^ BG(/, J; K) is a basic graph and Q is given by 

Q = {^i.i^fi\-nK) (64) 

where Tf^ is a right inverse of the block matrix Fj^j, that is ^iIj^^k ~ ^ik- 

Proof. In the appendix, we prove Lemma 1161 which states that for a graph F with input 

vertices /, output vertices J and measuring vetices K , the identity 

n[r\iK] = r(F|(?^'^)w(-F^j(?^'^|g'^)*U[p|p„c] . (65) 

holds if the equation 

r'K^jq'^'+p"' = (66) 

is valid. If we assume that F is basic, then for each measurement outcome p^^ the equation 
(|66|l has a solution q^'-' which is given by 

q""' = -ffflp'"" (67) 
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where Tfj^ is a right inverse for Fj^-j which exists since the block matrix is surjective 
(note that F]^^ is injectivc). By inserting this solution into 1)68(1 we obtain 

U[r|/K] = r(F| - f f//^')w(F^,f &"^| - f ,W^)*U[r|p.K] • (68) 

Thus we conclude that 6 = {^Kj^fKl ~ ^ik) indeed compensating. Vice versa, if we 
assume that there exists a compensating map, then U[r|/K] is an isometry. By Lemma 1101 
which is also proven in the appendix, it follows that F is basic □. 

As already mentioned, for a basic graph F G BG(/, J; K) the isometry \1[t\ik] is the encod- 
ing isometry of the graph code associated with F [Sl^]. Now we conclude from Theorem |51 
Corollary 6 Each stabilizer code can be implemented on a one-way quantum computer via 
an algorithm which is composed of the elementary steps )00 IT). 

4-3 Standard measurements and operations on graphs 

As long as we are concerned with basic graphs F G BG(/, J; K) the statement of Theorem|3 
tells us that the randomness of a measurement in x-basis at IK can be compensated via a 
conditional local unitary operation. In other words, we can solve the "first task" that goes 
along with one-way quantum computing. Thus for the implementation of the algorithm Tr 
it remains to tackle the second task which consists in computing the standard measurement 
operations. For this purpose, we simplify the channel Tr by removing the measurement 
vertices K such that we are left with input and output vertices only. In fact, it turns out that 
the standard measurement operations at K can be interpreted as operations on graphs. 

For a basic graph F G BG(/, J;if) we call a subset N C JK of measuring and output 
vertices pre-removable with respect to F if the block matrix is invertible, where denotes 
the inverse. We write PR(/, J, K\r) for the collection of all pre-removable sets. If a pre- 
removable set N contains no output vertices, we call it removable. 

Let M be a finite set which contains the same number of elements as J n A^. Let z^: M — > 
J n A^ be a bijective map. Then we denote by P the connecting graph on (J n N)M which 
connects the vertex m G M with z^(m) G J H A^ by one line. We obtain a new graph F + ? 
on the vertices IJKM which is, in particular, a basic graph with input vertices /, output 
vertices (J \ N)M and measuring vertices NK . For this graph the block matrix 

{T + dr^=V%+1)^=T^^ (69) 

is invertible. N is removable for F -f i?, and we can build the Schur complement 

X^(F + v)=T%- {T% + vf^^)f^^{V% + v-ir) (70) 

which is a graph with vertices NM , where N is the complement of A^ in UK. 

An important observation is here, that, if the underlying graph F is basic, then the graph 
Xjv(F + '!)) is also basic. Namely, for any basic graph A with input vertices /, output vertices 
J and measuring vertices K the following holds: Suppose that A^ C A" is removable. The 
Schur complement XjyA is a graph with input vertices /, output vertices J and measuring 
vertices M ~ K\N . In order to show that XjvA is a basic graph we study the kernel of the 
block matrix (XatA)/!/ which is nothing else but analysing the system of equations 

K\.JMTM_K.]MTM xJMJNiN „IM_n /71 \ 

[^nJ^)im<1 ~^im9 - ^n^iaiI -U (71) 
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which can equivalently be written as 

A/^^(9^^^-AM,g^^^) = 0. (72) 

Since A is basic graph, it foUows: If ql^ belongs to the kernel of A;J^, then Kjf^ = 
implies that = 0. Thus q'^' = follows from (EH if 9^ = q^*^- Aj^Af^^g^*^ is contained 
in the kernel of A^^. Indeed we find that 

Afx(<z"^ - Aj^Af,,g^^^) = Af,,/^^ - Aj^Aj^Af,,/^^ = (73) 

holds. Thus we have shown that the Schur complement of a basic graph is again a basic 
graph. Now we can apply this fact to the graph K = T which proves that XAr(r + P) is 
indeed basic. 

The prescription F i— > Xjv(r + 9) is an operation on graphs which leaves the number of 
input vertices as well as the number of output vertices fixed, because |Jn = \M\ and 
thus |J \ iV| + = I J| - IJA^I + \M\ = \ J\. On the other hand, the number of measuring 
vertices is reduced by |iCnA^|. We shall see, that standard measurement operations in x-basis, 
performed at a subset iYnif of measuring vertices K, can directly be viewed as applying the 
map F I— > Xat (F + v) to the underlying graph. 

Theorem 7 Let F G BG(/, J; K) he a basic graph, let N G PR(/, J, K\r) be a pre-removable 
set and let v. J C] N ^ M be a bijective map. Then the identity 

^Rfnivl "[r|/i^l = 1^ U[x„(r+9)|/K\jv] (74) 

holds for some phase k, where F^^m is the local Fourier transform with respect to the con- 
necting graph v acting at the output vertices J H N . 

Proof. We first apply Lemma Hlsl given in the appendix, to the basic graph T + V which 
has input vertices J, output vertices (J \ N)M and measuring vertices NK. The set N is 
removable for the graph F + P and the operator 

^ir+v\iNK]^[XN{r+v)\iK\N] e CI/ (75) 

is a multiple of the identity in 21(7) since it can be shown that the operator H75|l commutes 
with all Weyl operators with help of Lemma El Thus we conclude that there is a phase k 
such that 

^[r+v\iNK] = 1^ ^[XN{r+v)\iK\N] (76) 

holds. Note that k has to be a phase since both operators, u^r+viiNK] and U[Xjv(r+P)|/if\Af] 
are isometrics. Finally we employ Lemma ^1 which states that the identity 

F[r)M]U[Y\iK] = U[r+v\iKN] (77) 

holds. This means that the composition of the Fourier transform -F^m] with the isometry 
upi/x] can be viewed as joining the vertices of the connecting graph D to the vertices of 
F and adding the edges of both graphs. Combining the equations (|76|) and H77|) yields the 
desired identity ifTH) □. 

The identity (|74|) simplifies, if we just consider a removable set N e PR(/, J, K\^)- It is 
not necessary to make use of a bijective map ly and the Fourier transform can be skipped. 
The Schur complement for F is then given by: 
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Corollary 8 Let T e BG(/, J; K) be a basic graph. If the set N C K is removable, then the 
identity 

U[r\iK] = K U[x„r|/K\Ar] (78) 

holds for some phase k. 

ThcorcmClcan be used to remove all measurement vertices in a systematic way by just con- 
sidering one- and two-elementary sets. The following proposition is a further useful ingredient 
for this purpose. 

Proposition 9 Let T G BG(/, J, K^) be a basic graph. Then for each measuring vertex n € K , 
one of the following statements is true: 

(i) The vertex n has a self-link, V{n,n) ^ which implies that the set {n} is removable. 

ill) The vertex n has no self-link, T{n,n) = but it is connected to another measuring vertex 
n' e K, i.e. T{n,n') ^ 0. This implies that the set N = is removable. 

(iii) The vertex n has no self-link, T(n,n) = but it is connected to an output vertex j S J, 
i.e. T{n,j) ^ 0. This implies that the set N = is pre-removable. 

Proof. Wc choose a measuring vertex n e K. The statement corresponding to case (i) is 
obvious. Suppose now that m is a measuring or an output vertex in JK which is connected 
to n. This implies indeed that the set N = {n, m} is removable since the determinant of the 
block matrix is non-vanishing dct(r]^) ~ T{n,n)T{m,m) — r(ri, m)^ = — r(n, m)^ 7^ 
and it follows that is invertible, i.e. the statements (ii) and (iii) hold. To complete the 
proof we have to exclude the case where the measuring vertex n is connected to non of the 
measuring or an output vertices. In this case we find for a register configuration q^^ = 7^ 0, 
which only has a non-zero entry at position n, that the identity Tj^q^^ = V'l^q"' = holds. 
This contradicts the assumption that F is a basic graph, i.e. the kernel of Tj^ is trivial □. 

We can derive from the results of the previous discussion an inductive algorithm for re- 
moving the measuring vertices of a basic graph F £ BG(/, J]K). Suppose after s steps, we 
have removed the measuring vertices Mg such that we are left with Ks = K \ Mg remaining 
measuring vertices and a basic graph Fj € BG(/, Js, Kg). The set of output vertices Jg has 
the same cardinality as J and consists of a part Jg n J of initial output vertices and a part 
Js \ J of exchanged output vertices. Then we remove a further measuring vertex n g Kg 
conditioned by one of the possibilities, given by Proposition |5| 

We also have to keep in mind that some operations on graphs has to be compensated by 
a product of local Fourier transformations Fg-. L2{¥'^') L2{¥'^) which operate trivially on 
the qudits with positions J H Jg. 

(i) If the vertex n has a self-link, Tg(n,n) ^ 0, it is removable and we build the basic graph 

F,+i := X„F, (79) 

with input vertices /, output vertices Js+i = Jg and measuring vertices A'^+i :~ Kg \ n. 
The Fourier transform i^s+i = Fg remains unchanged. 

For the binary case, the Schur complement (|79|1 can be obtained by a useful graphical 
role: Take the subgraph by removing the vertex n. Add further edges to it by mutually 
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connecting those vertices which are hnked to n, including self links. Note that edges are 
added modulo 2. Fig. ^ represents this procedure. 

(ii) The vertex n has no self-link, but it is connected to another measuring vertex k € K^, 

i.e. rs(7i, k) 7^ 0. We build the basic graph 

Ts+i Xin^fcjTs (80) 

with input vertices /, output vertices Jg+i = Js and measuring vertices Kg^i := Kg \ 
{n,k}. The Fourier transform Fg^i = Fs remains unchanged. 

Again we have a graphical role for obtaining (|8()|l for the binary case. By following the 
strategy to remove all measurement vertices with self-links first, we may assume without 
loss of generality that both vertices n and k have no self-links. Then wc can proceed 
as follows: Build the subgraph by removing the vertices n and k. Then add edges by 
connecting each vertex, which is linked to n, with every vertex that is linked to k (See 
Fig. El for illustation) . 

(iii) The vertex n has no self-link, but it is connected to an output vertex k Cz Js- Let m be 

an additional vertex. Then 1™ + is the adjacency matrix of the graph that connects 
the vertex I with m by one line. We build the basic graph 

r,,+i :=X{„,,.}(r,, + i;" + i^) (81) 

with input vertices /, output vertices Js+i = (Js \ k)m and measuring vertices K^+i 
Kg \ n. The new Fourier transform is given by Fg+i — FgFll'^]. 

Concerning binary systems, we may use here the same graphical role as described for 
case (u). 



Example I: Fig. [7| shows a simple example for a qubit system (d = 2). The initial graph 
Fq has input vertices {1,2}, measuring vertices {3,4,5,6}, and output vertices {7,8}. The 
adjacency matrix is given by 
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(82) 



where the lines separate input, measuring and output vertices. The corresponding graph is 
shown in Fig. |3 Clearly, the block matrix 



{3,4,5,6|7,8} _ 
^{1,2|3,4,5,6} ~ 
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StepO 




Fig. 7. Step 0: The initial graph Fq. Step 1: The set measuring vertices {3,4} is removable. 
According to (ii), we obtain a new graph Fi with two connected measuring vertices {5, 6}. Step 2: 
The set measuring vertices {5,6} is removable. According to (ii), we finally obtain the graph F2 
which has no measuring vertices left. 



has maximal rank which imphcs that To is a basic graph. The measuring vertex 3 fulfills the 
condition (ii). It has no self-link but it is connected to the measuring vertex 4. Thus the set 
{3,4} is removable and we can build the Schur complement 
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(84) 



which is a graph on the remaining vertices {1, 2|5, 6|7, 8}. This operation is illustrated by 
Fig.d Step 1. 

Now, the measuring vertex 5 satisfies condition (ii) . It has no self-link but it is connected 
to the other remaining measuring vertex 5. The set {5, 6} is therefore removable with respect 
to the graph Fi and we can build the Schur complement as well 



X{5^6}ri = 
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0/ 



(85) 



which yields the final graph r2 that only lives on the input vertices {1,2} and the output 
vertices {7,8} (Fig.d Step 2). 



Example II: Fig |S1 shows a basic graph F which represents the underlying entanglement 
resource for running the algorithm under consideration. For illustration, we choose here an 
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algorithm which is encoded by measuring the qubit at position {1} in x-basis, the qubit 
locahzcd at {2} in z-basis, and the qubit sitting at {3} in y-basis. This scheme corresponds 
to a measurement graph A which is also depicted in Fig. |H1 




Fig. 8. Step 0: The left figure shows a graph F with measuring vertices {1,2,3} representing 
the underlying multipartite entanglement recource for performing algorithms. The right figure 
represents a measuring graph A which corresponds to the measurement strategy we wish to perform 
at the outputs {1,2,3}. {4} labels the auxiliary vertex and syndrome vertices are labeled by 
{5, 6, 7}. We perform local measurements in x-basis at {1}, in 2:-basis at {2} and in y-basis at 
{3}. This yields, combined with the corresponding measurement graph, the graph Fq. Step 1: 
The measuring vertices {2,4} of Fq are removable. According to (ii) we obtain the graph Fi. 



Step 2 




Fig. 9. Step 2: The measuring vertices {3} of Fi has a self-link and can be removed by (i) which 
yields the graph F2. Step 3: The measuring vertices {1} of F2 has a self-link and can also be 
removed by (i) which finally yields the graph F3. 

According to Theorem |31 we can substitute the graphs F and A by one graph Fq = F — 
Z{5 g 7}A which represents an equivalent entanglement recource on which only measurement 
operations in x-basis has to be performed (Fig. |H1 Step 0). First we remove, as shown by 
Fig.lHl Step 1, the vertices {2,4} by making use of (ii) which yields the graph Fi. This graph 
differs from the underlying graph F in the following manner: The measuring vertex {3} has 
an additional self- link and the vertex {2} is removed from the cluster. In view of that we can 
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formulate the following role: 

• A measuring vertex, which is measured in .x-basis, remains unchanged. 

• A measuring vertex, which is measured in z-basis, is removed from the cluster. 

• A measuring vertex, which is measured in y-basis, gets an additional self-link and is 
measured in x-basis. 

The measurement vertices {1,3} of the graph Fi can completely be removed by applying the 
strategy, which is described above. The vertex {3} has a self-link and can immediately by 
removed by (i) as shown by Fig.lHl Step 2. Again we obtain a graph r2 with one remaining 
measuring vertex {1} that has a self- link. Thus we can apply the strategy (i) once more which 
produces the final graph (Fig. Step 3) 

The present example shows one interesting feature which goes along with the removing 
of measuring vertices. The underlying graph F represents a configuration of qubits that are 
arranged in a cubic lattice with next-neighbor Ising interaction. Measuring qubits at positions 
{1, 2, 3} in X-, Z-, and y-basis respectively is equivalent to preparing a configuration of qubits, 
according to the graph F3, that evolve via this pattern of two-qubit Ising interactions and 
self-interactions. This can be interpreted as simulating the interaction, given by F3, by using 
next-neighbor Ising interactions and by performing local measurement operations. 

5 Conclusion 

We close the present paper by discussing an idea what the "essence" of an algorithm is. The 
procedure, given in Section |31 by a finite number of steps )00)01 nl), describes how an 
algorithm can be implemented by local measurement operations on cluster states. In the first 
step, the information is encoded into the cluster via preparing the input qudits at positions / in 
a state that carries the information we wish to process. The remaining qudits are individually 
prepared in the standard state. Within the next step we apply one elementary step of the 
dynamics which is given by a weighted graph F on the union of the set / of input vertices, the 
set O of output vertices, and mutually disjoint sets of measuring vertices {Mt\t = 0, • • • ,n}. 

In the first measurement round, we measure the input qudits in x-basis and we perform a 
local graph measurement on the qudits which are localized at M = Mq. The corresponding 
measurement graph A lives on the union of the set M = Mq, a set K of auxiliary vertices, and 
a set L of syndrome vertices. After this procedure has been performed, we obtain a random 
measurement result {p^,q^) (Fig. EJ. 

By making use of Theorem |21 only measurements in x-basis are needed in the first mea- 
surement round: We measure the qudits which are localized at the input vertices / and the 
measuring vertices KM in x-basis. This produces a random measurement outcome p^^^^ . 
Provided F — Z^A is a basic graph with input vertices / output vertices J and measuring 
vertices MK, Theorem |S1 can be applied which states that each random measurement out- 
come p^^'^^ can be compensated by applying a phase space translation that operates at the 
positions J. This translation corresponds to the vector (a'^|&'^) € 'B.'^ which is explicitly be 
given by 

= F^,,(r^ZlA);*V'' and = (F - Z^A);,,,,/^^^ . (86) 
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Here (F-ZlA) 

IKM ^ right inverse for the block matrix (F — Z^A) 



IKM 
JM ■ 



r-z,A 




^0 




^0- 




Fourier transform 



Fig. 10. Removing measuring vertices: Prepare the qudits at JKM individually in the standard 
state, apply an elementary step of the dynamic associated with V — Zx,A, measure the qudits at 
IKM in x-basis. Then perform, conditioned by the measurement outcome, a local unitary opera- 
tion at J that compensates the measurement result. This operation coincides with a operational 
scheme that does not make use of graph measurements: Prepare the qudits at Jf. individually 
in the standard state, apply an elementary step of the dynamic associated with V^, measure the 
input qudits at / in x-basis, finally perform, conditioned by the measurement outcome, a local 
unitary operation at J that compensates the measurement result. 



With help of Theorem [7| it fohows that the measuring vertices KM can be removed 
successively from the graph F — Z^A which we choose as initial graph Fq — F — Z^A with 
input vertices /, output vertices Jq — OMi ■ ■ ■ Mn and measuring vertices Kq ~ MK. After a 
finite number k of steps we obtain a graph F^ that just have input vertices /, output vertices 
Jfc and a local Fourier transform Fk that operates on the qudit positions Jfc \ ( Jo H Jfc). Thus 
the channel that is represented by Fig. can be realized by the simplified scheme which is 
represented by Fig. 1101 The local Fourier transform Fk can be absorbed into the compensating 
part and it is the "essence of the algorithm" what is left: The graph F^ with input vertices /, 
output vertices O and measuring vetices Mi, • • • , M„. However, the graph F^ is not unique 
and depends on in which order the measuring vertices are removed. 
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Appendix A 

Stabilizer codes and algebras 

In Section[21we discuss the general theory of stabilizer codes. From the mathematical per- 
spective the basic object is the abelian C*-algebra, the stabilizer algebra, 21(7, J, K\A), which 
is naturally represented on the Hilbcrt space L2{V'^). This part of the appendix serves several 
useful lemmas that enable to perform an explicit decomposition of the natural representation 
of 21(7, J, K\A) into irreducibles. 

The basic objects, which arc concerned here are the operators 

V[A|q^] = Vd''^^' <i>ljKL]Ui^)'^lJK]Tl[qL] (A.l) 

which play a central role for obtaining an explicit decomposition of the natural representation. 
In addition to that we are also concerned with the closely related operators 

U[r|p/K] =\/d'^^'$[p,K]w(r)$[Kj] (A.2) 

which describe the operation which has been performed after the measurement in x- basis 
produces the random result p^^ . The operator U[r|/i<-] = U[p|o//f] corresponds to the standard 
measurement outcome. It is an isometry whenever T S BG(/, J; K) is a basic graph, i.e. the 
block matrix Tj^ is injective: 

Lemma 10 // and only ifV is a basic graph in BG(/, J; K), then U[p|/^] is an isometry from 
L,{¥',) tntoL^iWi). 
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Proof. In order to decide whether U[r|7K] is an isometry, we now study the operator 
Ujp|j^,jU[r|/x]- Its action on a function -0 e L2{¥^j) is given according to 

(ufp|,,^]U[r|,K,V)(g2') = rf"'" / dqi Kiqi\qi)4>{qi) (A.3) 

where the "integral kernel" K is given by 

K(<72'k() = / dq^dq^dq' r{-r\qi^ + q-')T{T\qi^ + q-') (A.4) 

The map ty '■ q^ T{T\q) is a project representation of the additive group F^"^^ which satisfies 
the identity T{T\qi +92) = T{T\qi)T{T\q2)x{qi\^q2) ■ Making use of this fact, equation (|A.4p 
turns into the following form: 

Kfe'kf) = T{-T\qi)T{T\q{) I dqf dg2^r(-r|#)r(r|qf ) 

X X(rf 92^1 - q^)x{rfqi\q^) I dq'\{q' \Ti^{q{^ qi'<)) . (A.5) 



The integration over the variables q"^ can explicitly been performed where we keep in mind 
that the Fourier transform on a finite abelian group yields J dq'^x{p'^\q'^) ~ S{p'^). Thus we 
conclude that the integral kernel K fulfills 

Kfe'kO = ri-T\qi)rir\qi) J dgf d#r(-r|#)r(r|gf ) 

X Xirfqil q^)xiTfq{\q^)Sirj^{ql'' - qi'')) . (A.6) 

We exploit the invariance of the Haar measure dg^ and perform the variable transformation 
I3 = q2 — qi which yields 

K{qi\q{) = T{-^T\qi)T{T\q{) j dq^dq^ r{^T\qf )x{Tfqi\~qf) 

X x{Tf{q{-qi)-Tlq^\q'^)5{Ti{q{-qi)-VU?) (A.7) 

where we have used the identity T(-r|gf + )r(r|gf ) = T(-r|gf )x(rf | - ). Now 
we perform the integration over the variables q^ and we find that the equation 

Kfe'kO - T{-T\qi)T{T\q{) I dq^ r(-r|gf )x(rf ^2' | - q^)5{Ti^iqi - qi - q^)) (A.8) 

is valid. If we assume that F is a basic graph, then Tj^ (9i — ~ ^f^) = implies that q( = q^ 
and q^ = hold. Inserting this into IjA.Sp yields via performing the integration over q^ : 



n<l'M) = d-\^\5{qi-qi). (A.9) 
Finally, we insert this identity into Equation (|A.3|) and we obtain 

{ul^^jj,^U[rliK]i^){qi) - dl^l / dqi Siqi - qi)^{q{) = ^{qi) (A.IO) 
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which imphcs that u^ri/K] is an isomctry. On the other hand, one sees from equation (jA.Sp 
that the function K cannot be supported on the diagonal ql = (/I if there exists a non-zero 
vector q^^ with Tj^q^^ = 0. Hence, if F is not basic, then U[r|/K] cannot be an isomctry □. 

The following lemma presents a useful relation between the isometrics vjai^lj and the 
operators U[r|p/]. 

Lemma 11 Let T be a graph with input vertices I , output vertices JM and let K be a graph 
without input vertices, output vertices M and measuring vertices K . Then the identity 

^[A\q^]^[T\p'\ = ^[qL]U[T-A\QMK^pI^ (A. 11) 

holds for all register configurations q^,p^ . 

Proof. We insert the definition of the operators V[yv|gi] and U[p|p/] and we compute the 
expression 

^[Ak^]"[rb^l = ^/^'^'^'''^[,L]*fKM]^(-A)*[/^L]*[p^]"(^)*[./M] • (A.12) 

The operator V[A|gt] does not operate on the qudits that are localized at input vertices /. 
As a consequence, the co-isometry ^Jj,/] commutes with V[A.|gi]- Moreover, the operator u{T) 
operates trivially at the positions KL and commutes therefore with ^[kl]- Keeping in mind 
that |M| = |L|, this yields 

(A.13) 

where we have used the fact that u{T — A) = u(r)u(— A) holds which is a direct consequence 
of the definition □. 

Lemma 12 For all classical register configurations q^ , the linear operator v^^^^l-^ is an isom- 
ctry which maps the input Hilbert space L2(F^) into the output Hilbert space L2{¥'^). 

Proof. The operator Hj^t] is an isomctry and it is therefore sufficient to show that the 
operator 

= ^ <f*lIKL]^{miJK] (A.14) 

is isometric. We may interprete the vertices IL as "inputs" and K as "measuring" vertices. 
According to the assumption (G2), the block matrix A/^^ is invertible, in other words, 
A is a basic graph in BG{IL, J; K). Thus Lemma EH can be applied which proves that 
V[A|,i] = U[A|/KL]n[,i,] is an isometry □. 

Before the next lemma is presented, we note here a useful observation which will be used 
several times within this appendix. For any subset N C M the vector ^[^v] is invariant under 
all shifts x((7^) and ([n] is invariant under all multiplier z(p^). Moreover, according to the 
definition of x- and z- basis, the identities x{q^)([^-^ = C[qN] and z(p^)^[Ar] — ^[pW] are valid. 
This implies the useful relations 

w(/^k^On[A.] = n[,«]w(/AA^|/AA^) (A.15) 

Mp''\<l"')^[m = <i>[,«]w(/AA^|^M\A.) (A.16) 
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Lemma 13 The range of the isometry V[A.|gi] is the multiplicity space of the character w^qL] 
of the stabilizer algebra 2l(/, J, A'| A) which is uniquely determined by its the expectation values 
on Weyl operators according to 

W[,.]{MA'jK<l'''W')) = T(A|g^^^) , (A.17) 
where {Ajj^q"^^ Iq"^) is in the isotropic subspace S[jj^fc\A]- 

Proof. For proving the eigenvalue equation (|24|l . we consider a register configuration q'^^ 
which belongs to the kernel of the block matrix Aj^ . With help of (jA.lSp and HA.16|I we 
verify that the identity 

c V[A|,^] = $fj^^i,u(A)x(g'^^^)$[j^]n[i] (A.18) 

holds, where c is an appropriate normalization constant. In the second step, we make use 
of the Weyl commutation relations as well as the fact that u(A) implements the symplectic 
transformation {p'^^\q^') {p^' - Aq*^*^ with M = IJKL. This gives 

c V[A|,.] = r(-A|g^^)w(A;;^,g-^^|g'^)$*_^,^^,,_g,j,.(A)z(A^^g-^^) 

X x(/)$[j,^]n[i] . (A.19) 

We exchange the order of z(Aj^(7'^^) and x((7^), making use of the assumption that q"^^ 
belongs to the kernel of Aj^ and applying the relation (|A.15|I once more which yields the 
identity: 

c V[M,.] = r(-A|g'^^)x(-A^^g'^^|/)w(A^^./^|g^)$f,^^jH(A)$[,;^]n[,.j . (A.20) 

According to the definition of the isometry V[A|qt], we find 

V[M,.] = r(-A|g'^^)x(-A5^g'^^|/)w(Ai^./^|g'^^)v[A|,.] . (A.21) 

Since there are no edges between syndrome vertices we conclude T{A\q'^^)x{Ajj^q'^^ \q^) 

coincides already with T{A\q'^^^) which proves the lemma □. 

Lemma 14 The family {'V[\^qL]\q^ E F^} of isometrics is mutually orthogonal 

^[Mgf 1^[A|,,-] = HQi - qi) 1/ (A.22) 

and complete 

^[Ak^lV[A|,^] = 1./ ■ (A.23) 

Proof. We compare the expectation values of two characters '^[p^] and n7[pi,j for a given Weyl 
operator w(Aj^g'^^|(7"') which belongs to the stabilizer algebra (see HA.lTfl ) by analyzing the 
ratio of both eigenvalues 

t i\ = "-lT$^xiI^i''iqi 9.^)k'^^) = x(Ai^fe^ - ■ (A.24) 

T{A\q'^^+q^) r{A\q^) 

Here we have used the fact that r(A|g|^) = T{A\q^) = 1 since there are no edges which 
connect syndrome vertices. The fraction of the eigenvalues ljA.24p is = 1 for all configurations 
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q.iK ^ kcrAj^ only if Ajf^ {q^ — g^^) lies in the range of A/^, i.e. there is a configuration q^^ 
such that 

Ai^(9f - q^) + A/f = (A.25) 

holds. By assumption (G2), the matrix Aj^j^ is invcrtible which implies q^^ = and q^ = q2 ■ 
Thus the characters ^[q^\ and ^[q!^\ coincide only if q^ = 92 ■ Lemma [T51 the ranges of the 
isometries vj^Yj^i,] and Vj^j^i,] are therefore orthogonal to each other if q^ and g|' are distinct. 

The range of each isometry V[A|qt] has dimension rf'^L Wc have already shown the or- 
thogonality ljA.22p which implies that the projection, defined by the left-hand side of (|A.23p . 
projects onto a rfl^l+l-^l dimensional space. By assumption (Gl) we have |/| -I- \L\ = | J| and 
(f03|) follows □. 

Appendix B 

Operations on cluster states 

In this part of the appendix we develop the tools for tackle local measurement operations 
in x-basis (this implicitly covers also measurements in y- and z- bases) which are applied to 
cluster states in an appropriate manner. A Weyl operator, which acts on the input Hilbert 
space ^2(1^^) commutes with the isometry \1[t\ik] according to the following useful relation: 
Lemma 15 Let T be a graph with input vertices I , output vertices J and measuring vertices 
K . Then for each register configuration p^ the commutation relation 

t{TW"^) M^r\iK] Mp'W) = ^{VijKl'^'W') ^[T\IK] (B.l) 
holds if q'^^ is a solves the equations 

v'^AjW" (B.2) 

0^ = ^lu<^'"' . (B.3) 

Proof. We first recall the fact that the isometry ^'yjK\ is invariant under all shift operators 
x(g'^^) (see (IA.16|) ) and that the unitary operator M(r) implements the symplectic transfor- 
mation (j>^\q^'^^) ^ {p^ — Tq^'^^ Iq^"^^), i.e. the corresponding commutation relation 1)12(1 
holds. From this we conclude that 

T{T\q"'<) u{T)^[jK] w(/|g^) = w(/ - Tq'''^\q'''^^)uir)<^>[jK] (B.4) 

is valid. Now we apply the co-isometry '^^ both sides of the Equation l|B.4p which yields 

r{r\q"'')u[riiK] w(/|g^) = w(-r,Vg^'^^|g-^)u[p|r.K^,..^_p.] . (B.5) 

Here we have used the identity ((A.16|l . Therefore the commutation relation (|B.1() follows if 
P' - ^ijki"^ = is fulfilled □. 

A similar result is given by the next lemma which can be used to prove the existence of 
compensating maps (see Subsection 14. 21 for the definition of this notion). 

Lemma 16 Let T be a graph with input vertices L , output vertices J and measuring vetices 
K , then the identity 

U[r|/K] = r(r|g^-Ow(-ri.,;g^'^|g-^)*U[p|p.K] . (B.6) 
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holds if the equation 

(B.7) 

is fulfilled. 

Proof. Due to the shift invariance of the isometry ^[kj], i-c. I)A.16(I . the identity 

dl^^lupip.K] - <I>fp„,]u(r)x(g^^)<I>[,„] (B.8) 

is fulfiUed for all classical register configurations g^''. We exchange the order of the operators 
u{T) and x(g^'^) which yields by keeping (|A.16p once more in mind: 

dl^^lu[r|p.K] = r(-r|g^''^)w(-r^^g^'^|q-')<I>fp„,+p^.^^,,,]u(r)$[K,/] . (B.9) 

This implies that (|B.6|I is true if the classical register configurations p^^ and q^'-' solve the 
equation T^j^jq^'' + p^^ = □. 

A heuristic analysis of standard measurement operations in .T-basis suggests that they are 
directly related to opeartions on graphs. The following two lemmas serves as key ingredients 
for proving that this intuition is indeed correct. 

Lemma 17 Let T e BG(/, J; K) be a basic graph and let N G K be removable. Then the 
identity 

Tq^'^^ = (Xjvr)g^'^*^ (B.IO) 
holds, M = K \ N , for all q"^^ which solve the equation 

r\^Kq"''-p' (B.ll) 
for some p^ ,q^ G F^. In particular, it follows that 



iT\q''^)=T{XNT\q'''') (B.12) 



is true for all q^'^^ that fulfill (|B.11|) . 

Proof. We first split the system of equations (|B.lip into two parts 



r|)W^ + rjv^ V = / (B.13) 

TU"" + ^^jM<l'''' =^ (B.14) 



which is equivalent to 



T\%q"''+T^^''q''=p' (B.15) 

^f^T^jMq"'^ (B.16) 

since, by making use of the assumption that N is removable, has an inverse f j^. Thus we 
can substitute q^ into the first equation which yields 

p\7A/ TJM __T^IM TJM j^IMf^Nj^N JJM _T m 1 7^ 



Finally one verifies that (|B.10p holds by making use of (|B.16p via inserting it into the left hand 
side of IjB.lOp . Equation ljB.12|) follows immediately from the definition of the phases T(r|-) 
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Lemma 18 Let T e BG(/, J; K) be a basic graph and let N C K be removable. Then the 
operator 

u*r|/x]U[x„r|/x\Ar] G CI7 (B.18) 
is a multiple of the identity in 2l(/). 

Proof. By Lcmnia ITHl the isometrics U[-p\iK] and U[x„r|/j\f]5 M ~ K\ N, commute with 
a Weyl operator 'w{p^\q^) according to l|B.ip . Since F is a basic graph, the block matrix 
Tj^ is surjective and for each p^ there is q^'^^ which fulfills the equations p^ — Tjq^ = 
T'jj^q-^^ and T^kI'^^ = This implies that the equations p' = {XNTYjjj^jq^-"^' and 

{yi.M^)YjMl^'^^^ ~ are valid due to the fact that is invertible. Hence we conclude that 
the relation 

T(r|<7^'^^) U[r|,;^] w(/|q^) = ^{-TijKq'"'\q') U[r|/K] (B.19) 

as well as 

r(~X^r|<z^'^*^) ^{p'Wr ufx„r|/M] = U[x„r|/M] O^nTYj j Y (B.20) 

is fulfilled. Due to Lemma [T7I we also conclude that the identities T{y^M^W^'^^) ~ '''(rl'?'^''^) 
and \>f{{y^M^Yi.jMq''^^W^) = '^{^ijkI^'^^Ii'^) are also true. Combining both equations 
(|RT9)| and (|R^ implies that 

U[x„r|/A/]U[r|/x] w(p^|90 = w(/|g^) ufx„r|/A/]U[r|/K] (B.21) 

holds, i.e. ^'^XNr\iM]^U^\ii<] commutes with all Weyl operators. Thus, it is indeed a multiple 
of the identity □. 

A further important relation shows that the application of a Fourier transform to a cluster 
state can also be interpreted as an operation on the underlying graph by adding additional 
vertices and edges. 

Lemma 19 Let F G BG{I, JN; K) be a basic graph with output vertices NJ, J D N = (d. 
Moreover let v: N i—t M be a bijective map and v the connecting graph that connects each 
vertex n ^ N with its counterpart v{n) G A/. Then T + v is a basic graph in BG(/, JM; KN) 
and the identity 

^[p^-r]U[r|/A'] = ^[T+V\IKN] (B.22) 

holds, where F[i}^'] is the local Fourier transform associated with the connecting graph v. 

Proof. According to our assumption F is a basic graph to which we add the connecting graph 
V which yields the graph T+v with input vertices /, output vertices JM and measuring vertices 
KN . This is indeed a basic graph since the kernel of the block matrix {J^ + vYink^ — ^Ynk'^ 
is trivial. This can be observed as follows: Suppose that (F + vYifix^ q^^^ ~ holds. 
Then this implies q^ = and therefore q^ = 0. Thus we conclude that ^^'lYkq^^^ — 
Y\jfqiK = is valid. But since F is basic, we find = 0. 

Applying the isometry U[p|/^] to a vector iIj G L2(]F^) can explicitly be computed according 

to 

(U[r|/x]^)(9''^) = c J dg^^r(F|g"^^)V'(g^) (B.23) 
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where the constant c = \/d ' ' ' ' is for normaHzation. Concatenating this operator with 
the Fourier transform i^rpj/i for the connecting graph 9 yields 



= cc'J dq'^^ x(/^|^>i^^9^)r(r|g™)^(g^) (B.24) 

with c' = \fd ' '. Making use of the definition of the bi-character x('l') ^^^d the phases t(-|-) 
we verify that the identity 



x(/^|Pjyg^)r(r|g--«) = exp 



IJNK\ 



a 



exp 



cxp 



T{T + V\q) (B.25) 



holds for all register configurations q = qU^^t^ . By inserting this relation into the Equation 



(|B.26ll and by using the definition of the isometry VL[r+v\iKN] gives rise to the desired identity 



which completes the proof of the lemma □. 



(B.26) 



